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SECTION  I 


r 


INTRODUCTION 

Under  certain  rotation  rates,  a  rotor,  supported  by  an  oil-lubricated 
journal  bearing,  can  experience  excessive  shaft  vibration  (Fig.  1).  Two 
types  of  oscillations  can  exist  called  oil  whirl  and  oil  whip  (Ref  1-6). 

Oil  whirl  produces  a  circular  orbit  in  a  direction  of  the  same  sense  as  the 
shaft  rotation  but  at  a  value  of  about  45%  of  the  shaft  rotational  speed 
(Fig.  2).  Oil  whip  is  a  more  violent  oscillation  which  occurs  at  the  reso¬ 
nance  frequency  of  the  shaft  system.  Although  this  phenomenon  has  been 
investigated  for  some  time  and  identified  as  a  self-excited  oscillation  due 
to  viscous  effects,  no  satisfactory  analytic  prediction  method  has  been 
developed.  It  is  the  purpose  of  this  note  to  provide  such  an  analysis. 
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SECTION  II 


GOVERNING  DYNAMIC  EQUATIONS 


Consider  an  unbalanced  rotor  with  mass,  m,  shaft  elasticity,  k, 
rotating  at  speed,  to,  under  the  influence  of  fluid  forces,  X  and  Y 
(Ref.  1,  2). 


Hence 

2 

mx  +  kx  =  mn<o  cos  co  t  +  X 

2 

my  +  ky=mno)  sin  ui  t  -  mg  +  Y 
The  classical  method  (Ref.  2)  for  describing  the  damping  forces  is 


X  =  -cx 
Y  =  -cy 


The  two  dynamic  equations  are  uncoupled  and  can  be  readily  solved  as  follows 

c  2  k 

(for  small  damping,  in  which  (— )  <<  — )  : 

Im  m 


x 


ct 

2m 


cos 


2 

(w  t  +  <(>)  +  r  cos  uit 

o  Z  Z 

OJ  -U) 

o 


=  r  e 
o 


ct 

2m 


sin  (w  t  + 
o 


<to  + 


mo 


2  2 
(0  -to 

o 


sin  wt 


where 

2  k 

u>  =  — ;  r  and  6  are  ascertained  from  initial  conditions, 
o  mo 

These  results,  however,  give  little  insight  into  the  cause  of  oil 
whirl  or  whip.  The  reason  for  the  failure  is  believed  to  be  the  inability 
to  adequately  describe  the  fluid  forces.  Therefore  another  mathematical 
model  for  the  prediction  of  these  forces  was  developed. 
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SECTION  III 


GOVERNING  FLUID  MECHANIC  EQUATIONS 


In  this  section  the  fluid  motion  between  two  non-concentric  rotating 
cylinders  will  be  analyzed  to  obtain  the  fluid  dynamic  force  on  the  inner 
cylinder  representative  of  a  journal  bearing  (Fig.  3).  The  radius  of  the 
inner  cylinder  is  'a'  and  the  radius  of  the  outer  cylinder  is  ’a+b1, 
where  ^  <<  1.  To  analyze  this  problem  two  coordinate  systems  must  be 
defined,  i.e.,  an  inertial  system  which  is  fixed  at  the  center  of  the 
outer  cylinder  (x,y)  and  a  rotating  coordinate  system  with  a  moving 
origin  referenced  to  the  center  of  the  inner  cylinder  (x ’ ,  y').  The  set 
of  polar  coordinates  associated  with  these  two  systems  are  (r,0)  and  (r',6') 
respectively.  The  inner  cylinder  is  considered  to  be  rotating  at  a  con¬ 


stant  value  of  O  about  the  center  of  the  outer  cylinder  with  an  orbit 
radius  of  eb.  The  incompressible  Navier-Stokes  equations  in  polar  coor¬ 
dinates  for  this  rotating  system  are  as  follows: 

3u  ,3  /  t  \  n 
— ,  +  (r  v)  *  0 


Du  uv 
Dt  r' 


3p  ,  /r72  u  23v 

jpf5T  +  v(7u-~-j--5— 

r'  r  39 


)  -  29v 


Dv  u  3p  . 

dF  •  T7  =  -  +  V(v 


2  v 
v  - 


?  ? 

r'  r'  30' 


)  +  0 (2u  +  r'0) 


where 


D _ 3_  3_  3 

Dt  3t  v  3r '  U  r ' 30 ' 


+  i  1_  + 


8r’2  r’  3r'  r'230'2 


Consistent  with  lubrication  theory  these  equations  may  be  greatly  simpli¬ 
fied  for  small  gap  thicknesses. 


Let  r' 


a+z 


where 

z 

a 

<<  i 

and 

_v 

u 

«  i 

where  z  has  limits  between  zero  and  h,  the  gap  thickness.  Using  the  law 

of  cosines  a  relationship  for  h  may  be  obtained. 

(a+b)^  =  (be)^  +  (a+h)^  -  2be  (a+h)  cos  (u-6') 

For  the  case  at  hand,  —  <<  1,  hence 

a 

h(8')=b  (1-e  cos  6')  is  a  good  approximation. 

An  order  of  magnitude  analysis  of  the  Navier-Stokes  equations  produces 
the  classical  lubrication  equations. 


Q  = 


$ 


udz  =  flow  rate 


dp  3  u  ,. 

__r_  -  p  — -  =  pressure  gradient 

3z“ 

with  boundary  conditions 
u(6' ,  z=0)  =  (w-0)  a 
u(9 ' ,  z=h)  =  -  0a 

The  pressure  must  satisfy  a  periodic  boundary  condition. 

p(e')  =  p(e'+2ir) 


Integrating  the  pressure  gradient  equation  twice  with  respect  to  z 
produces  a  relationship  for  u. 

u  =  f  +  C  z  +  C, 

aud6  212 

Applying  the  two  boundary  conditions  for  u  on  the  cylinder  walls 
eliminates  the  two  integration  constants. 


u+8a 


I^dF  z(z'h)  -  T  (z“h) 


A 


Note  that  from  this  expression  the  shear  stress  on  the  wall  of 
the  inner  cylinder  can  be  obtained  as  follows: 

3u  ,  h  dp  pauj 

T=p  —  (z=0)  =  ~  2 ^  -  IT 

Integrating  u  with  respect  to  z  produces  the  flow  rate. 


Q  = 


h 

s 


udz  =  <f  -  B)  ah  -  ~ ^ 


or  rearranging  this  equation. 


bdp/d9’ 


1-a-e  cos0’ 


12a2p(^  -  0)  (1-e  cos  0')^ 


where 


a  = 


ab(|  -  6) 

This  equation  can  be  integrated  with  respect  to  0'  to  determine  the  pressure 
distribution  around  the  cylinder. 


Ap  = 


O’ 

5 


d£ 

de 


,  de* 


The  following  integral  is  required:  (Ref.  7) 
0' 


5 


(l-ct-e  cos  9')  ,Q,  -e  a  sin  9*  e(6-a)  sin  6 ' 

(1-e  cos  0’)  2 (1-e  ) (1-e  cos  0’)  2 (1-e  )  (1-e  cos  0*) 


_  0  * 

(e-e2  a)  -1  ,/l-e2  tan  2  x 

+  - 9  c'Jj  tan  ( - r~ - ) 

(1-e2)572  1-6 

where  8=2  (1-a-e2) 

The  periodic  boundary  condition  on  pressure,  i.e.,  p(0')  ■  p(9?  +  2x) , 

requires  that  the  coefficient  of  the  last  term  must  vanish.  Therefore, 

2 

8-e  a=0,  provides  a  relationship  between  the  flow  rate  (Q)  and  rotational 

speed  (0) , 


5 


1 


2(l-e  )  _  Q 

2  +  e2  ab(|  -  0) 

Using  this  condition  to  eliminate  a  and  8  produces  the  final  relation- 
ship  for  the  pressure  distribution  over  the  cylinder  (called  the 


Harrison  Equation,  Ref.  8-12). 
Ap 


12  a2(|  -  9)e  sin  6'  (2_e  cos  6') 


2  2 
b  (2+e  ) 


(l-e  cos  8  ’ ) 2 


One  additional  integration  will  produce  the  desired  relationship  for  the 
forces  on  the  inner  cylinder. 


F  =  -aL 
r 


2it  2n 

Ap  cos  8'  d8'  -aL  ^  t  sin  0'  d0'  =  0 
o  o 


$ 


-aL 


2  tt 

S 


Ap  sin  0'  d0'  -aL 


2n 


5' 


cos  0'  d0'  = 


24irua  Le 

2  2  t - 

b  (2-H;  )  /l-e 


2  l<«-  f>  +|d-3a)(0-  f)  +  Jof] 


where  o(e)  *  — 2+e? ^  (1  -  /l-e2) 
6e 


For  -  «  1 
a 


24iruaJLe 


b2(2+e2)  /Hi 


2 1*  -!  (i  -  °i>] 


The  small  term,  o—  in  the  above  expression  for  Fa  is  the  only  remaining 
a  o 

contribution  of  the  shear  stress  integral  indicating  that  the  pressure 
term  dominates. 
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Having  obtained  a  description  of  the  fluid  forces  on  non-concentric 
rotating  cylinders,  we  can  now  return  to  the  dynamic  equations  for  the 
rotor  system. 
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SECTION  IV 


DYNAMIC  EQUATIONS  IN  POLAR  FORM 


Since  we  found  it  convenient  to  express  the  fluid  forces  in  polar 
co-ordinates,  we  choose  to  similarly  convert  the  dynamic  equations. 
Utilize  the  following  transformation. 


Let  x  = 

r  cos 

0 

*  Fr 

=  X 

cos  9 

+ 

Y 

sin 

0 

y  = 

r  sin 

0 

5  F0 

=  Y 

cos  0 

+ 

X 

sin 

0 

The  governing  dynamic  equations  in  polar  form  become 

"2  2 

m(r-r  0  )  +  kr  *  mnw  cos  (0-wt)  -  mg  sin  0  +  F^ 

•  •  ••  2 

m(2  r  0  +  r  0)=  -mnu  sin  (0-ut)  -  mg  cos  0  +  F 

V 

These  equations  are  non-linear  and  must  be  linearized  to  continue  further. 
Consider  only  small  perturbations  in  both  r  and  0 
r  =  be  +  6 (t) 

0  *  ft  +  q(t) 

From  the  analysis  of  the  fluid  forces  we  determined  the  following 

F  =  0 
r 

F.  =  m  SI  be  C  (0-ft) 

0  q 

u  n  u)  bx  _  w 

where  ft  =  -r  (1-  o— )  _  -r 

l  a  2 


and  C 

q 


3 

•247Tga  L 


mftb3(2+e2)  /l^e 


2 


Note,  the  effect  of  the  shear  stress  term  (<*~)  in  the  expression 
for  ft  reduces  the  value  of  ft  to  slightly  below  50%  of  u>.  This  result  is 
in  qualitative  agreement  with  the  experimental  data  (Ref.  2). 

Using  these  relationships  and  neglecting  higher  order  terms  produces 
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the  following  set  of  linear  equations. 

fi  +  (--  -  ft^)  6-2ftbeq  =  noi^  cos  (n-u)t-g  sin  ft  t+F  /m-(—  -  ft^)be 
m  r  in 


2ft6  +  beq  -  (ftbe  C^)q  =  -  nu  sin  (ft-w)t-g  cos  £2  t 
The  right  hand  sides  of  these  equations  contain  the  forcing  functions  and 
lead  to  the  particular  solution.  The  homogeneous  solution  contains  the 
source  of  the  self-excited  oscillation  and  can  be  obtained  as  follows: 

Let  =  AeXt 
be 


and  ^  =  Be 


Xt 


Setting  the  right  hand  sides  to  zero  we  obtain 

-2X 


2  k  2 

xz+(-  -  ftZ) 

m 


X-ftC 


X 


2ft  q 

This  produces  a  cubic  characteristic  equation. 

3  2  k  2  k  2 

X  -ft  C  Xz  +  (-  +  3ft  )X  -  ft  C  (-  -  ST)  *  0 

q  m  q  in 

Routh's  discriminant  states  that  no  positive  real  root  exists  provided 
all  coefficients  and  determinants  of  this  equation  are  positive.  Hence, 
stable  motion  occurs  provided 

1.  C  <  0 

q 

k  2 

2.  -  -  ft  >  0 

m 
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SECTION  V 


ANALYSIS  OF  RESULTS 


In  the  previous  section  two  conditions  were  derived  for  stable  motion 
of  the  journal  bearing.  If  the  two  conditions  are  satisfied  then  a 
steady  autorotation  will  persists  at  0=0=0. 5  w  and  r=be,  producing  a 
circular  orbit.  This  motion  is  descriptive  of  oil  whirl  for  which  experi¬ 
mental  evidence  (Ref.  2)  indicates  a  stable  autorotation  at  about  0'  0.45  w. 
Our  predictions  indicate  that  C  is  always  negative. 


C 

q 


247rpa3L 

mOb3(2+e2) 


<  0 


Hence,  the  first  condition  is  always  satisfied.  Stable  rotation  is  a  con¬ 
sequence  of  the  fact  that  dF./d8  C  is  negative,  forcing  the  rotation  rate 

0  q 

to  return  to  the  autorotation  value  (Q  =  |)  if  perturbations  occur  on  either 

side  of  this  value.  Also  equilibrium  is  attained  at  zero  torque  or  F  =  0 

which  occurs  when  6  =  ^  .  Also,  experiments  show  that  oil  whirl  occurs 

only  if  fi  <  /k  (Fig.  2)  which  verifies  the  second  condition, 
tn 

If  the  second  condition  is  violated  i.e.,  Q  >  fk  ,  structural  diver- 

m 

gence  of  the  shaft  occurs  in  which  the  centrifugal  force  of  the  whirl  becomes 
greater  than  the  elastic  restoring  force  of  the  shaft.  The  eigenvalue,  X, 
becomes  positive  indicating  a  departure  solution.  The  journal  deflection 
will  become  excessive  and  strike  the  bearing  wall  during  rotation  creating 
unusual  journal  position  patterns.  This  motion  at  the  divergence  state  is 
descriptive  of  oil  whip. 
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SECTION  VI 


CONCLUSIONS 

Self-excited  oscillations  in  journal  bearings  called  oil  whirl  and 
oil  whip  have  been  predicted  by  obtaining  an  analytic  description  of  the 
fluid  forces  and  inserting  these  into  the  dynamic  equations  for  the  rotor 
system.  Oil  whirl  is  shown  to  be  an  autorotation,  f2,  of  the  journal 
centroid  at  a  rate  of  one  half  the  shaft  rotational  speed.  This  value  is 
in  good  agreement  with  experiment.  Stable  rotation  is  a  consequence  of 

the  fact  that  dF./d0  is  negative  forcing  the  rotation  rate  to  return  to 

y 

the  autorotation  value  (^)  if  perturbations  occur  on  either  side  of  this 
value . 

Oil  whip  is  shown  to  occur  when  the  frequency  ">f  autorotation  exceeds 

the  natural  frequency  of  the  shaft  system,  i.e.  fi  /k  . 

m 

At  this  value  structural  divergence  occurs  and  the  centrifugal  force 
exceeds  the  elastic  restoring  force.  Positive  eigenvalues  indicate  de¬ 
parture  solutions  which  are  in  agreement  with  experiment. 

Although  simplifying  approximations  to  the  Navier-Stokes  equations 
produced  satisfactory  results,  it  is  recommended  that  additional  research 
be  accomplished  to  assess  the  consequences  of  this  mathematical  model. 
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